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ABSTRACT 

/  /1^‘LO 

Heat  transfer  and  pressure  data^ obtained  with  sharp  flat  plate  models  at 
zero  and  large  angles  of  attack  in  the  CAL  hypersonic  shock  tunnel.ar-re— p-re-sented-. 
These  high  Mach  number  data  to  24)  extend  from  the  classical  thin 

boundary  layer  regime  to  near-free -molecule  conditions  and  are  discussed  within 
the  framework  of  existing  theory.  The  large  angle-of-attack  results  verify  viscouf 
shock-layer  theory  and  define  the  low-density  conditions  where  transport  effects  at 
the  shock  wave  _tsh-Ock-wave  slip)  first  become  important. 

The  data  obtained  at  zero  angle  of  attack  are  compared  with  theory  to  define 
the  fluid  mechanism  that  governs  the  low-density  effects.  The  strong  shock-wave 
approximations,  the  vorticity  interaction,  and  the  shock-wave  heating  effect  are 
aliown  lu>  be  unimportant  in  this  case.  It.  is~sJh  own -that  for  the  present  experiments 
the  effects  of  transport  processes  at  the  shock  wave  are  numerically  small  in  com¬ 
parison  with  surface  slip  for  the  sharp  flat  plate  at  zero  angle  of  attack,  and  it  is 
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^postulated  that  surface  slip  is  the  dominant  low-density  effect. 

The  slip  boundary  conditions  are  discussed  to  demonstrate  that  there  is  a 
pressure  jump  at  the  surface,  analogous  to  the  classical  velocity  and  temperature 
jump.  This  effect  is  evaluated  using  existing  theories  and  is  in  qualitative  agree¬ 
ment  with  the  experimental  data.  The  slip  boundary  conditions  are  examined  in 
light  of  the  present  experiments  at  zero  angle  of  attack,  and  it  is  shown  that  the 
usual  first-order  boundary  conditions  are  inapplicable  for  the  case  of  a  cold  wall_ 
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Constants  determined  by  the  Prandtl  number  (Eq.  (1)  ) 
Free  constant  in  the  velocity  jump  and  enthalpy  jump 
Chapman -Rube  sin  constant 
Modified  Chapman-Rubesin  constant, 

Specific  heat  at  constant  pressure 
Stanton  number 
Enthalpy 

Thermal  conductivity 
Eq.  (30) 

Mach  number 
Pressure 

Reservoir  pressure 
Pitot  pressure 
Prandtl  number 

Rate  of  heat  transfer  per  unit  area 
Gas  constant, 

Reynolds  num 
Temperature 

Modified  reference  temperature  (Eq.  (5)  ) 

Velocity 

Mean  molecular  velocity 
Distance  measured  from  plate  leading  edge 
Distance  measured  normal  from  plate  surface 
Angle  of  attack 
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Eq.  (3) 

Ratio  of  specific  heats 
Eq.  (46) 

Boundary-layer  displacement  thickness  without  surface  slip 

Perturbation  boundary-layer  displacement  thickness  due  to  slip 

Density  ratio  across  a  shockwave 

Molecular  mean-free-path 

Viscosity 

Density 

Shock  wave  angle  measured  from  surface 


Eq.  (5) 
Eq.  (27) 
7  Eq.  (12) 


Subscripts 

C  Equilibrium  nozzle  expansion 

n  Nonequilibrium  nozzle  expansion 

O  Reservoir  conditions 

I  Conditions  at  the  base  of  the  continuum  boundary  layer 

W  Conditions  at  the  wall 

00  Conditions  in  the  free  stream 

3  Conditions  behind  shock 
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I.  INTRODUCTION 


The  recent  interest  in  manned  re-entry  vehicles  has  precipitated  a  con¬ 
certed  attempt  to  describe  the  transition  between  continuum  fluid  dynamics 
and  free-molecule  mechanics.  The  problem  arises  from  the  need  during 
manned  re-entry  to  maintain  the  deceleration  forces  tolerably  small  while  the 
high  kinetic  energy  associated  with  orbital  or  space  flight  is  dissipated  in  the 
atmosphere.  This  dictates  that  much  of  the  atmospheric  flight  take  place  at 
high  altitudes  where  the  ambient  density  is  low,  in  a  transitional  regime  that 
falls  between  classical  continuum  flows  and  free-molecule  flows. 

One  aspect  of  this  transitional  flow  regime  is  its  complexity.  If  one 
adopts  either  the  continuum  or  the  kinetic  viewpoint,  it  is  necessary  to  con¬ 
sider  a  variety  of  phenomena  introduced  or  influenced  by  the  details  of  the  con¬ 
figuration.  The  sharp  flat  plate,  which  is  employed  in  the  present  research,  is 
an  ideal  research  tool  because  the  governing  low-density  phenomena  can  be 
isolated  and  studied  with  relative  ease.  Competing  low-density  effects  intro  ' 
duced  by  the  configuration  details  (for  example,  shock-wave  curvature  effects 
due  to  nose  bluntness)  are  suppressed  with  a  flat  plate. 

The  problem  of  a  sharp  flat  plate  in  low-density  supersonic  flow  has  been 

considered  over  the  past  decade.  The  earliest  theoretical  work  was  that  of 

1,2  3  4  5 

Shen,  followed  by  Lees  and  Probstein,  Lees,  Li  and  Nagamatsu, 

Stewartson,  Kuo,  and  Cheng.®’  ^  The  problem  has  been  experimentally  in¬ 
vestigated  at  low  Mach  numbers  by  Schaaf,  et  al^  and  by  Bertram,**  and  high 

w  ,  12 
Mach  number  investigations  have  been  made  by  Hammitt  and  Bogdonoff  and 

9,  1  3 

by  Hall  and  Golian.  As  a  result  of  these  and  other  investigations,  the 

viscous  interaction  problem  on  a  flat  plate  is  well  understood  in  both  the  weak 
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and  the  strong  interaction  regimes  as  long  as  the  Mach  number  and  Reynolds 
number  are  large  enough  to  suppress  the  low-density  effects  neglected  in  the 
viscous  interaction  regimes. 

These  additional  low-density  effects  cause  substantial  departures  from 
the  thin  boundary  layer  results  as  the  density  is  decreased.  The  departures 
have  been  qualitatively  anticipated  on  the  grounds  that  as  the  leading  edge  is 
approached,  the  free-molecule  limit  must  be  obtained.  However,  there  has 
been  uncertainty  in  defining  the  fluid  mechanisms  that  would  govern  the  initial 
departures,  in  specifying  a  valid  theoretical  model  for  use  in  describing  the 
flow  field,  and  in  determining  the  magnitude  of  these  departures. 

The  purpose  of  this  paper  is  to  present  experimental  heat  transfer  and 
pressure  data,  obtained  with  sharp  flat-plate  models  in  the  CAL  hypersonic 
shock  tunnels,  which  extend  from  the  thin  boundary  layer  limit  to  near-free- 
molecule  conditions.  The  intent  is  to  demonstrate  the  magnitude  of  the  low- 
density  departures  from  thin  boundary  layer  theory,  to  test  the  available  low- 
density  theories,  and  to  identify  the  dominant  fluid  mechanism  causing  the  initial 
departures  from  the  thin  boundary  layer  regime. 

In  this  connection,  the  applicable  theories  are  reviewed  in  Section  II  to 
establish  a  basis  for  interpreting  the  experimental  results.  The  classical  slip 
boundary  conditions  are  reviewed  there  to  show  that  there  is  a  pressure  jump 
at  the  surface  which  is  larger  than  the  usual  slip  effect.  Sections  III  and  IV 
describe  the  experimental  apparatus  and  test  conditions,  and  the  experimental 
results  are  presented  and  discussed  in  Section  V.  The  first-order  slip  boundary 
conditions  are  appraised  in  Section  V  to  show  they  are  inapplicable  for  the  case 
of  a  cold  wall. 


2 


The  results  for  the  flat  plate  at  large  angles  of  attack  are  discussed 
first.  This  class  of  flows  is  of  particular  interest  since  only  one  low-density 
mechanism  (the  transport  processes  at  the  shock  wave)  is  important.  Thus  it 
is  possible  to  evaluate  the  magnitude  of  this  mechanism  and  to  show  that  it  is  a 
small  effect  at  zero  angle  of  attack.  The  more  complicated  case  of  a  flat 
plate  at  zero  angle  of  attack  is  discussed  next.  The  indications  are  that  the 
dominant  low-density  mechanism  in  this  case  is  slip  at  the  surface. 
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II.  STATUS  OF  THEORIES 


There  is  a  large  body  of  theoretical  literature  dealing  with  the  sharp  flat 
plate  in  the  thin  boundary  layer  regime.  A  comprehensive  review  of  the  prob¬ 
lem  can  be  obtained  from  the  literature  cited  in  the  Introduction.  In  addition,  a 

substantial  body  of  literature  is  developing  on  low-density  phenomena.  A 

15 

general  review  of  the  field  has  been  given  recently  by  Sherman.  The  purpose 
in  this  Section  is  to  set  out  the  theoretical  developments  pertinent  to  the  inter¬ 
pretation  of  the  experimental  results.  In  this  connection,  the  theory  for  wedge 

flows  will  be  reviewed  first  in  Sub-Section  A,  emphasizing  Cheng's  viscous 
1  7 

shock-layer  theory.  The  reason  for  this  choice  is  that  this  theory  provides  a 
solution  which  yields  the  Blasius  result  at  high  densities  and  approaches  the 
free -molecule  limit  for  vanishing  Reynolds  number. 

The  theoretical  results  pertinent  to  the  sharp  flat  plate  at  zero  angle  of 
attack  are  reviewed  in  Sub -Section  B,  beginning  with  the  thin  boundary  layer 
(weak  and  strong  interaction)  results.  One  approach  to  treating  low-density 
effects  is  to  include  additional  low-density  mechanisms  in  the  thin  boundary 
layer  results  with  a  perturbation  calculation.  Each  of  these  mechanisms  is  re¬ 
viewed  to  provide  a  basis  for  evaluating  their  effect.  The  mechanisms  consid¬ 
ered  here  are  1)  breakdown  of  the  strong  shock-wave  approximations,  2)  the 
vorticity  interaction,  3)  the  shock-wave  heating  effect,  4)  molecular  processes 
at  the  surface  (surface  slip),  and  5)  transport  processes  at  the  shock  wave 
(shock-wave  slip).  It  is  noted  that  effects  stemming  from  finite  shock-wave 
thickness  should  be  considered.  There  is  presently  no  theoretical  treatment  for 
the  sharp  flat  plate  known  to  the  authors  for  assessing  this  effect. 

Each  of  the  above  mechanisms  has  been  discussed  in  the  literature.  The 
vorticity  interaction  and  shock-wave  heating  effect  arise  from  the  assumption 
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that  the  local  inviscid  flow  properties  are  fixed  by  the  local  shock-wave  angle. 
Actually,  the  flow  is  processed  by  a  stronger  portion  of  the  shock  wave  near 
the  leading  edge,  and  the  local  flow  is  nonuniform  with  a  higher  temperature. 
The  effects  of  transport  processes  at  the  shock  wave  has  received  recent  atten¬ 
tion  and  become  important  when  the  shock  layer  becomes  fully  viscous.  Under 
these  conditions,  the  classical  (uniform  flow)  Rankine-Hugoniot  conditions  are 
inapplicable,  and  it  is  necessary  to  account  for  the  gradients  behind  the  shock 
wave. 

In  the  discussion  of  surface  slip,  some  detail  is  given  to  deriving  the  slip 
boundary  conditions  using  the  classical  Maxwell  model.  This  is  done  in  order 
to  show  that  there  is  a  pressure  jump  at  the  surface,  analogous  to  the  velocity 
jump.  In  addition  it  is  noted  that  the  classical  temperature  jump  is  inconsistent 
for  hypersonic  flows,  and  that  conceptual  consistency  requires  that  an  energy 
or  enthalpy  jump  be  considered. 

Another  approach  to  the  leading-edge  problem  on  a  sharp  flat  plate  is  to 
postulate  a  flow  model  independent  of  the  thin  boundary  layer  models,  and  to 
seek  a  solution  which  applies  at  the  leading  edge  and  approaches  the  thin  bound¬ 
ary  layer  solutions  in  the  high -density  limit.  The  theories  derived  from  both 
the  continuum  viewpoint  and  the  molecular  viewpoint  are  reviewed  in  the  latter 
paragraph  of  Sub-Section  B. 

(A)  Large  Angle  Wedge  Flows 

The  wedge  flow  problem  at  high  Reynolds  numbers  is  classical  and  has 
been  treated  with  the  assumption  that  the  rate  of  boundary  layer  growth  is  small 
in  comparison  with  the  wedge  angle.  In  this  scheme,  the  wedge  pressure  is 
determined  by  the  wedge  angle,  and  the  boundary  layer  characteristics  can  be 
calculated  using  the  Blasius  solution  based  on  conditions  behind  the  shock 
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waves.  ^  As  the  density  is  decreased,  one  might  anticipate  the  boundary  layer 


displacement  effects  would  become  important  for  certain  wedge  angles 

/  M2/1  1 

^  ^ |  J  .  These  first-order  effects  of  boundary  layer  displacement  have 

^  8  9 

been  theoretically  investigated  by  Cheng  ’  and  experimentally  verified  by  Hall 
9  13 

and  Golian  ’  as  a  part  of  a  more  comprehensive  study  of  the  combined  effects 
of  boundary  layer  displacement,  nose  bluntness,  and  angle  of  attack.  The  re¬ 
sults  approach  the  strong  interaction  result  in  the  limit  of  vanishing  wedge  angle, 
and  the  wedge  result  in  the  limit  of  high  density. 

The  low-density  effects  in  wedge  flows  have  been  extensively  investigated 

17 

in  a  study  by  Cheng  of  the  blunt- body  problem  at  low  Reynolds  number.  An 
important  feature  in  this  analysis  is  that  if  the  wedge  angle  is  sufficiently  large, 
one  can  suppress  the  effects  of  the  strong  shock-wave  approximations  and  the 
effects  of  shock-wave  heating  and  the  vorticity  interaction.  Further,  it  is  shown 
that  the  effects  of  surface  slip  and  those  associated  with  the  shock-wave  struc¬ 
ture  are  of  higher  order  than  those  due  to  transport  processes  at  the  shock  wave. 
This  class  of  flows  is  of  considerable  importance  since  it  isolates  a  single  fluid 
dynamic  mechanism.  By  examining  these  flows  in  theoryand  experiment,  it  is 
possible  to  determine  the  ambient  stream  conditions  where  the  transport  effects 
become  important.  The  experimental  data  reported  in  Section  V  will  serve  to 
verify  the  theory  and  to  demonstrate  the  magnitude  of  the  effects  of  transport 
processes  at  the  shock  wave. 

(B)  Zero  Angle  of  Attack 

The  sharp  flat  plate  at  zero  angle  of  attack  in  low-density,  hypersonic 

flow  has  been  investigated  for  over  a  decade,  and  theoretical  solutions  are  pres- 

4  9 

ently  available  for  treating  the  boundary  layer  displacement  effects.  ’  A 
comprehensive  review  of  this  problem  and  the  theoretical  treatments  may  be 
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found  in  these  papers  as  well  as  in  Ref.  14.  The  approach  generally  employed 

P 

is  to  consider  two  regimes;  one  of  weak  interactions  (~  0  { 1 )  )  and  one  of 

p  Ro 

strong  interactions  (~  >>  1}.  The  subsequent  discussion  will  retain  this 

rao 

subdivision,  but  it  should  be  noted  that  local  similarity  applies  in  both  regimes, 

and  that  the  subdivision  is  necessary  only  because  of  the  approximations  used 

to  describe  the  pressure  ratio  across  the  shock  wave. 

Weak  Interaction  Regime.  There  are  several  solutions  available  for  the 

weak  interaction  regime.  The  one  used  here  is  that  given  by  Hayes  and  Prob- 
14 

stein.  Following  this  development,  the  surface  pressure  is,  to  the  first- 
order  in  the  viscous  interaction  parameter. 


P 

Poo 
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X  =  M: 


f, 


(1) 


where  C  is  the  Chapman- Rubes  in  constant  relating  viscosity  to  temperature. 
For  a  Prandtl  number  of  0.  725,  A  =  0.  968,  and  B  =  0.  145.  In  calculating 
the  heat  transfer,  it  is  found  in  Ref.  14  that  there  is  no  first-order  effect  due 
to  the  self-induced  pressure  gradient  and  that,  to  the  first  order  in  jC  >  the 
heat  transfer  is  identical  to  the  zero-pressure  gradient  result.^  Similarly, 
the  effect  of  the  induced  pressures  on  the  boundary  la>er  is  to  cause  an  over-all 
decrease  in  the  thickness,  but  there  is  no  first-order  effect  on  the  boundary 
layer  slope.  Consequently  it  is  concluded  that  the  weak  interaction  description 
of  the  pressure  distribution  is  accurate  to  the  second  order  in  X  ■  anc*  t'lat 
this  result  is  valid  for  X  £  JSiY-n 

^A(Pr)^eCp.) 

It  is  interesting  to  note  that  the  contribution  of  the  second-order  term  in 
the  pressure  is  quite  small  at  the  theoretical  limit  of  validity  of  weak  interaction 
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theory;  about  15%.  Consequently,  only  the  first-order  result  for  the  pressure 
distribution  will  be  used  here.  In  addition,  it  should  be  noted  that  if  the 
Howarth-Dorodnitsyn  transformations  are  applied  to  the  governing  boundary 
layer  equations,  it  can  be  demonstrated  that  local  flat-plate  similarity  applies, 
and  that  the  error  associated  with  neglecting  the  pressure  gradient  terms  is  of 

order,  £  —  Hence,  the  heat  transfer  rate  in  the  weak  interaction 

8  +  1 

regime  can  be  calculated,  assuming  a  Prandtl  number  of  one,  using  the  general 

9 

relation  for  heat  transfer  given  by  Cheng. 


P 


where 


I3 


2  L  Ti 


ACRO+2 


The  relation  in  Eq.  (3)  subsequently  will  be  compared  with  experimental  re¬ 
sults. 

As  noted  earlier,  there  are  a  number  of  solutions  available  for  the  strong 

interaction  regime.  The  solution  used  in  the  present  investigation  is  that  given 
9 

by  Cheng.  The  advantages  in  this  solution  are  in  the  fact  that  it  clearly  reveals 
the  role  of  wall  temperature  in  the  strong  interaction  regime  and  that  the  solu¬ 
tion  takes  account  of  the  effects  of  shock-wave  curvature  on  shock-wave  strength 
by  using  the  Busemann  formula  for  the  pressure  ratio  across  the  shock  wave. 
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(4) 


The  results  for  the  surface  pressure  and  heat  transfer  are 


£  (o.gg4+  i.73 ^-)m3G  =  aaLiiffi'fltJ 
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where  — 
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_  r  _  X-l 

Rb».  >  tf-H 


It  will  be  noted  that  the  Chapman- Rube  sin  constant  in  Eqs.  (4)  and  (5)  is  based 
on  a  mean  temperature  different  from  that  normally  used.  That  aspect  of  the 
problem  has  been  critically  examined,  ’  and  it  was  concluded  that  this  defi¬ 
nition  yielded  a  better  approximation  to  the  average  temperature  in  the  boundary 
layer. 

More  recently,  solutions  have  been  generated  for  the  boundary  layer  die- 

18  19  20 

placement  effects  using  high-speed  computers.  ’  y’  In  Ref.  18,  the  local 
similarity  concept  is  applied  throughout  the  weak  and  strong  interaction  regimes 
and  the  limits  of  validity  are  established.  The  complete  tangent -wedge  for¬ 
mula  (as  opposed  to  the  limiting  values  often  applied  in  weak  and  strong  inter¬ 
action  theory)  for  the  pressure  ratio  across  an  oblique  shock  wave  is  applied  in 
the  inviscid  region,  and  numerical  solutions  are  presented  for  a  variety  of 

cases.  The  theory  has  been  tested  against  existing  theories  and  experimental 

1 8 

data  and  found  to  be  in  good  agreement. 

The  Strong  Shock-Wave  Approximations.  The  accepted  practice  in  strong 
interaction  developments  is  to  assume  that  the  shock  wave  is  strong.  A  factor 


in  the  lower  limit  of  validity  on  ^  is  in  these  approximations.  Some  insight 
into  these  limitations  can  be  obtained  by  writing  down  the  oblique  shock-wave 
relations  for  the  pressure  and  density  ratio  across  the  shock  wave,  assuming 
the  wave  angle,  ,  is  small. 


The  strong  shock  approximations  and  their  limits  of  validity  are 
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for  >>  ('y-iKuO  (9) 


The  limits  of  validity  are  cast  into  this  form  since  the  pressure  is  measured 
in  the  experiments  and  consequently  provides  a  convenient  method  for  assessing 
the  strong  shock-wave  approximations.  It  can  be  seen  that  the  strong  shock- 
wave  approximation,  as  applied  to  the  pressure  across  the  shock  wave,  is 

P 

readily  satisfied  even  in  the  weak  interaction  regime  for  -pr  ">  1  However, 

roc 

this  approximation,  when  applied  to  the  density  ratio,  introduces  a  rather 
severe  limitation  which  suggests  that  the  strong  interaction  theory  should  not 
apply  for  pressure  ratios  of  less  than  about  25.  This  point  will  be  considered 
further  in  connection  with  the  experimental  results. 
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The  Vorticity  Interaction  and  Shock  Heating  Effects.  One  approximation 

in  strong  interaction  theory  is  to  assume  that  the  pressure  and  temperature  just 

outside  the  boundary  layer  are  equal  to  those  quantities  just  behind  the  shock 

wave  at  the  same  chordwise  station.  Actually,  the  streamline  which  intercepts 

the  edge  of  the  boundary  layer  at  this  station  passed  through  the  shock  wave  at 

a  more  forward  station,  and  hence  was  processed  by  a  stronger  portion  of  the 

bow  wave;  the  shock-wave  heating  effect. 

Another  assumption  made  is  that  the  velocity  is  uniform  in  the  inviscid 

flow  between  the  boundary  layer  and  shock  wave.  In  the  physical  situation,  the 

shock  wave  is  curved  near  the  leading  edge,  thereby  producing  a  nonuniform 

inviscid  flow  between  the  boundary  layer  and  the  shock  wave.  This  vorticity 

21 

interaction,  first  noted  by  Ferri  and  Libby,  is  important  in  that  it  alters  the 

boundary  condition  to  be  applied  at  the  edge  of  the  boundary  layer. 

The  effects  of  shock-wave  heating  and  the  vorticity  interaction  in  the 

22 

strong  interaction  regime  have  been  quantitatively  examined  by  Lees  and  sub- 
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sequently  refined  by  Oguchi.  In  Ref.  22  it  is  assumed  that  the  pressure  dis¬ 
tribution  is  given  by  strong  interaction  theory,  and  the  alterations  stemming 
from  these  higher-order  effects  are  examined  with  a  perturbation  calculation. 
One  assumption  made  is  that  the  changes  in  the  pressure  distribution  are  small 
in  comparison  with  the  strong  interaction  result,  It  is  demonstrated  that  the 
vorticity  at  the  edge  of  the  boundary  layer  13  invariant  with  chor 
while  the  shock-wave  heating  effect  is  governed  by  the  paramete 
The  quantitative  effects  on  surface  pressure  are  estimated  by  considering  the 
changes  stemming  from  the  alterations  in  the  boundary  conditions  at  the  edge  of 
the  boundary  layer.  It  is  demonstrated  that  for  an  insulated  plate,  both  the 
vorticity  interaction  and  the  shock-wave  heating  effect  are  governed  by  the 
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parameter 


and  that  the  effect  is  to  increase  the  pressure. 
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The  important  point  to  be  noted  here  is  the  form  of  the  result.  That  is,  within 
the  framework  of  this  perturbation  scheme,  the  influence  of  the  vorticity  inter¬ 
action  and  the  shock-wave  heating  can  be  examined  by  separating  the  strong 


interaction  effects  and  considering 
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as  influenced  by  the  parameter, 


Transport  Processes  at  the  Shock  Wave.  One  effect  which  could  account 

for  the  low-density  departures  from  strong  interaction  theory  is  the  effect  of 

the  transport  processes  neglected  in  the  strong  interaction  theoretical  model. 
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This  mechanism  was  originally  examined  by  Sedov,  was  subsequently  inves- 
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tigated  by  Probstein  and  Kemp  and  by  Cheng  in  connection  with  blunt-body 
flows  and  has  recently  been  included  in  a  formulation  of  the  sharp  flat-plate 
problem  by  Probstein  and  Pan. ^  The  problem  arises  in  writing  down  the  clas¬ 
sical  Rankine-Hugoniot  equations  for  the  shock  wave,  which  apply  if  the  shock 

layer  is  inviscid.  If  the  shock  layer  is  fully  viscous,  the  classical  relations  do 
14,  1  7 

not  apply.  '  However,  the  Rankine-Hugoniot  equations  can  be  generalized 
for  a  viscous  shock  layer  to  obtain  the  boundary  conditions  to  be  satisfied  at  the 
shock  wave.  These  can  be  cast  into  the  following  form  (Pr  =  1)  analogous  to 
the  surface  slip  boundary  conditions. 
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where  the  subscript,  S  ,  refers  to  conditions  behind  the  shockwave.  Using 
the  following  transformation, 

$r  -Sc,  ©  -  .  t  ■  /Wh >  T 

and  introducing  the  strong  interaction  result  for  the  pressure  (Eq.  (4)  ,  the 
boundary  conditions  at  the  shock  wave  (Pr  =  1)  are 
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It  can  be  seen  that  under  these  conditions,  the  parameter  governing  the  trans- 

fdp 

port  processes  at  the  shock  wave  is  |V|  \J  • 

Velocity  Jump  at  the  Wall.  Slip  at  the  surface  is  encountered  as  the  den¬ 
sity  at  the  wall  is  decreased  until  the  molecular  mean-free -path  at  the  wall 
becomes  appreciable  In  this  case,  one  must  take  account  of  the  molecular 
energy  and  momentum  exchange  at  the  wall  in  applying  a  continuum  boundary 
layer  analysis.  That  is,  the  velocity  at  the  base  of  the  boundary  layer  is  non¬ 
zero,  and  the  temperature  at  the  base  of  the  boundary  layer  can  be  different 
from  the  wall  temperature. 

The  classical  approach  to  the  velocity  jump  at  the  wall,  as  developed  by 
Maxwell,  effectively  is  to  regard  the  boundary  layer  as  a  nonuniform  continuum 
flow  except  for  the  layer  with  a  thickness  of  one  mean-f ree-path  at  the  wall.  The 
molecules  at  the  base  of  the  boundary  layer  are  assumed  to  have  a  Maxwellian 
velocity  distribution,  and  it  is  assumed  that  the  velocity  jump  is  given  by 
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where  X  is  the  mean-free -path,  bu.  is  a  constant  determining  the  magnitude 

of  the  velocity  jump,  and  the  subscript,  1  ,  refers  to  conditions  at  the  base  of 

the  continuum  boundary  layer.  The  constant,  bu  ,  is  determined  by  intro- 

P  >s 

ducing  the  tangential  momentum  accommodation  coefficient,  TV  >  and  con¬ 
serving  the  tangential  momentum  across  the  free-molecule  layer.  The 
conservation  of  tangential  momentum  is  expressed  as  follows: 


(15) 


The  term  on  the  left  is  the  net  tangential  momentum  transported  through  the 

boundary  layer,  the  first  term  on  the  right  is  the  tangential  momentum  transfer 

for  a  Maxwellian  gas,  and  the  second  term  on  the  right  is  the  additional  tangen- 
27 

tial  momentum  flux.  Substituting  Eq.  (14)  in  (15),  the  free  constant  is  ob¬ 
tained. 


(16) 


Enthalpy  Jump  at  the  Wall.  The  classical  statement  for  the  temperature 
jump  is  obtained  with  a  similar  computation  which  is  based  on  conserving  energy 
across  the  free-molecule  layer.  It  is  important  to  note  that  in  this  computation 
it  is  implicitly  assumed  that  the  gas  is  quiescent,  and  hence  the  computation 
does  not  rigorously  apply  for  the  hypersonic  case.  The  appropriate  statement 


v  The  tangential  accommodation  coefficient  is  defined  as 
is  the  tangential  momentum,  and  the  subscripts  t  and  f 
and  reflected  quantities . 
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refer  to  the  incident 
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is  that  there  is  a  jump  in  energy  or  total  enthalpy  at  the  wall  which  can  be  cal¬ 
culated  in  a  manner  analogous  to  that  for  the  velocity  jump. 


Hi  =  Hw  +  be  Xi  (  £  ^  )i 


(17) 


Introducing  the  thermal  accommodation  coefficient,  CLe  •  the  conservation  of 
energy  across  the  f ree-molecule  layer  is 


k(18) 


where  the  term  on  the  left  is  the  energy  transported  through  the  continuum 
boundary  layer,  the  first  term  on  the  right  is  energy  transferred  by  a  Maxwellian 
gas,  and  the  second  term  on  the  right  is  the  additional  energy  transported  to 
the  Maxwellian  gas.  Substituting  Eq.  (17)  in  Eq,  (18),  the  free  constant,  be  , 
is 


The  first  term  in  the  above  relation  is  that  obtained  in  the  classical  cal¬ 
culation  of  the  temperature  jump,  and  the  second  term  is  the  velocity  correction 
term.  The  usual  argument  for  neglecting  the  second  term  is  that  it  is  of  second 
order  in  the  mean-free-path  at  the  wall.  It  is  important  to  note  that  the  second 
term  can  be  neglected  regardless  of  the  magnitude  of  the  mean-free-path 

V- 1 

because  of  the  factor,  J - L  ,  Consequently,  the  enthalpy  jump  is  approximately 

0+  I 

equal  to  the  classical  temperature  jump  and  for  Pr  -  1  ,  1.0  ,  the  vel¬ 

ocity  jump  and  enthalpy  jump  are  equal.  The  heat  transfer  to  the  wall  is  the 


Defined  as  where  El  is  the  energy,  and  the  subscript,  W 

refers  to  the  wall  conditions. 
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term  on  the  left  in  Eq.  (18). 


This  statemen  of  heat  transfer  to  the  wall  is  comparable  with  that  given  by 
28 

Maslen,  taking  account  of  the  work  due  to  sliding  friction. 

Pressure  Jump  at  the  Wall.  The  velocity  and  enthalpy  jumps  are  the 

hypersonic  equivalents  of  the  classical  slip  boundary  conditions  usually  cited 
27 

in  the  literature.  It  is  important  to  note,  however,  that  the  energy  and  only 
one  component  of  the  momentum  have  been  conserved  in  deriving  these  rela¬ 
tions;  the  normal  component  of  momentum  has  not  been  considered.  In  writing 
the  conservation  of  normal  momentum  across  the  free-molecule  layer,  it  can 

be  shown  that  there  is  also  a  pressure  jump  at  the  surface.  The  surface  pres- 

.  14 
sure  is 
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(21) 


where  is  the  normal  momentum  flux  per  unit  area  incident  on  the  surface  and 
is  the  normal  momentum  flux  per  unit  area  of  the  molecules  re-emitted 


from  the  surface,  assuming  they  are  emitted  diffusely,  and 


■fit  i 


is  the  normal 


accommodation  coefficient.  In  the  spirit  of  the  classical  slip  model,  the  inci¬ 
dent  momentum  is  simply  that  of  a  Maxwellian  gas  with  an  additional  term  to 
account  for  momentum  transport  through  the  boundary  layer,  which  to  the  second 
order  in  is 
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R  = 


no  n, 
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(22) 


Defined  as  "fit  =  -2 - 

6  -  & 
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where  m  is  the  particle  mass,  LTT  is  a  mean  velocity,  and  H  is  the  number 
of  incident  molecules  per  unit  volume.  The  continuity  requirement  across  the 
free-molecule  layer  is 


m  Hw  iiw 


2X1  if\ ) 


(23) 


Combining  Eqs.  (21),  (22),  and  (23)  and  using  the  relations  in  Ref.  14  for  the 
Maxwellian  mass  flux  and  momentum  flux  it  follows  that  the  wall  pressure  is 


Within  the  framework  of  thin  boundary  layer  analysis,  it  is  usually  assumed 
that  the  normal  pressure  gradient  is  zero.  Equation  (24)  shows  that  in  this 
case,  there  is  a  pressure  jump  at  the  surface  resulting  from  the  temperature 
difference  across  the  free-molecule  layer.  The  importance  in  the  pressure 
jump  is  that  it  is  a  direct  effect.  The  velocity  and  enthalpy  jump  are  important 
in  that  they  alter  the  boundary  layer  displacement  thickness  which  in  turn  alters 
the  induced  pressure,  R  .  Because  of  the  pressure  jump,  the  wall  pressure 
is  different  from  the  induced  pressure  in  direct  proportion  to  the  velocity  and 
enthalpy  jumps.  Evidently  this  effect  has  not  been  previously  incorporated  in 
boundary  layer  analyses  that  include  surface  slip. 

In  concluding  this  discussion  of  the  slip  boundary  conditions,  it  should  be 
noted  that  the  classical  computation  reproduced  above  is  conceptually  incon¬ 
sistent  and  consequently  somewhat  unsatisfying.  In  particular,  in  writing  down 
the  energy,  mass,  and  momentum  balance,  it  is  first  assumed  that  the  gas  one 
mean-free -path  from  the  wall  has  a  Maxwellian  velocity  distribution,  implying 
that  the  boundary  layer  properties  are  constant  over  some  undefined  height.  The 
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fact  that  the  boundary  layer  is  non-Maxwellian  is  taken  into  account  by  adding 

on  additional  terms  proportional  to  the  continuum  gradients.  The  slip  boundary 

conditions  have  been  derived  by  Patterson  from  the  molecular  viewpoint  using 

29 

a  perturbation  scheme.  That  calculation  assumes  the  gas  is  slightly  nonuni¬ 
form  and  the  results,  given  to  the  first  order  in  \|  ,  are  in  close  agreement 

with  those  cited  here. 

Surface  Slip  in  the  Strong  Interaction  Regime.  The  classical  tempt... ature 
and  velocity  jumps  have  been  applied  by  a  number  of  investigators  in  analyses 
of  the  laminar  boundary  layer  at  low  densities.  The  general  approach  in  these 
calculations  is  to  determine  the  changes  in  the  boundary  layer  displacement 
thickness  due  to  slip,  and  the  attendant  changes  in  induced  pressure,  skin  fric¬ 
tion,  and  heat  transfer.  One  analysis  that  bears  directly  on  the  present  results 

30 

is  a  perturbation  calculation  by  Aroesty  based  upon  strong  interaction  theory. 
The  aim  there  is  to  examine  the  importance  of  surface  slip  in  the  strong  inter¬ 
action  regime.  The  analysis  is  based  on  the  classical  first-order  slip  boundary 
conditions  in  which  the  slip  velocity  and  temperature  jump  are  evaluated  at  the 
wall  conditions  as  opposed  to  conditions  at  the  base  of  the  boundary  layer,  i.  e. 


The  analysis  does  not  include  the  pressure  jump.  The  solution  is  carried  to  the 
point  of  setting  out  the  momentum  equation  describing  the  first-order  slip 
effects,  which  is  then  solved  approximately  for  the  case  of  an  adiabatic  wall. 
This  approximate  solution  is  based  on  the  assumption  that  the  primary  effect  of 
slip  stems  from  the  change  in  the  boundary  conditions.  The  result  obtained  for 
the  pressure  is 
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where  to  is  the  no -slip  strong  interaction  description  of  the  surface  pressure 

r* 

(for  example,  Eq.  (4)),  flo  is  the  no-slip  boundary  layer  displacement  thick¬ 
ness,  and  cO  is  the  change  in  displacement  thickness  stemming  from  the 
slip  boundary  conditions.  The  parameter,  u)  ,  in  its  general  form  is 


o)  = 


Xww  u  <30 

Z  £  u«  d*. 


where  Aw  is  the  molecular  mean-free -path  at  the  wall. 

By  invoking  the  strong  interaction  solution  for  the  pressure,  within  the 

>!< 

context  of  the  perturbation  analysis,  the  governing  parameter  becomes 

^  (2») 

For  an  adiabatic  wall,  it  is  demonstrated  that  the  effect  of  slip  is  to  decrease 
the  displacement  thickness,  and  the  final  result  for  the  pressure  is 

-§~  =  Po  I  ]  (29) 


The  analysis  concludes  that  for  a  Prandtl  number  of  one,  there  is  no  slip  effect 
on  the  heat  transfer  coefficient. 

Leading-Edge  Solutions.  A  number  of  investigators  have  examined  the 
low-density  leading-edge  problem  outside  the  framework  of  strong  and  weak 


There  appears  to  be  a  numerical  error  in  Aroesty’s  original  evaluation  of 
this  parameter. 
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boundary  layer  interactions.  The  approach  used  by  Oguchi  is  to  postulate  a 

wedge -like  continuum  flow  model  with  a  viscous  flow  filling  the  entire  region 

between  the  body  and  the  shock  wave.  A  solution,  assuming  slip  effects  at  the 

surface  to  be  negligible,  showed  that  a  limited  region  of  constant  pressure 

should  exist  at  the  leading  edge.  This  problem  has  been  re-examined  by  the 

same  author  using  the  classical  (velocity  and  temperature  jump)  first-order 

slip  boundary  conditions.  The  inclusion  of  the  classical  slip  boundary  conditions 

evidently  removes  an  anomaly  in  heat  transfer  that  appeared  in  the  earlier 

treatment,  and  the  heat  transfer  now  appears  to  approach  a  maximum  value. 

In  addition,  the  governing  parameter  that  arises  in  this  treatment  is  M 

Experimentally,  this  suggests  a  data  correlation  with  this  parameter  in  the 

strong  interaction  regime  and  beyond. 

The  leading -edge  problem  has  also  been  treated  with  kinetic  theory  by 
33  34  35 

Charwat  '  ’  seeking  a  solution  which  consistently  yields  the  f ree-molecule 

limit  and  defines  the  initial  portions  of  the  transition  to  the  continuum  regime. 

The  problem  is  complicated  since  within  the  framework  of  the  molecular  approach, 

it  is  necessary  to  define  the  velocity  distribution  function  for  the  molecules  after 

33 

each  collision.  The  uriginal  analysis  introduces  certain  simplifying  assump¬ 
tions  to  make  the  prjblem  tractable.  The  more  critical  assumptions  are: 
a)  the  emission  from  the  wall  is  equal  to  the  cross  stream  mass  flux  in  the 
free  stream  and  invariant  with  chordwise  iocation,  b)  the  gas  is  two-dimensional, 
and  c)  only  the  flux  of  fiee-stream  molecules  and  molecules  that  have  exper¬ 
ienced  one  collision  is  considered.  The  solution  for  the  surface  pressure  is 
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If  one  assumes  a  Maxwellian  distribution  of  velocities, 


The  heat  transfer  has  been  calculated  by  the  present  authors  from  the  original 
analysis  as 


In  a  subsequent  report,  the  original  solutions  were  modified  and  free 
constants  introduced  in  order  to  minimize  the  inaccuracies  caused  by  the  neces¬ 
sary  approximations.  These  constants  were  determined  by  assuming  a  normal 

momentum  accommodation  coefficient  of  unity  and  matching  the  solution  with 

3  6 

limited  experimental  pressure  data.  This  procedure  is  a  useful  device  for 

empirically  correlating  various  data.  However,  a  consistent  application  of  the 

determined  constants  in  a  relation  to  calculate  heat  transfer  leads  to  the  con- 

35 

elusion  that  the  thermal  accommodation  coefficient  must  be  0.  4.  Such  a  re¬ 
sult  is  contrary  to  the  experimental  data  cited  in  Section  V. 
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III.  EXPERIMENTAL  APPARATUS 


The  low-density  investigations  described  here  were  made  in  the  CAL  Six- 
Foot  Hypersonic  Shock  Tunnel.  Certain  high-density,  low  Mach  number 
experiments  were  made  with  the  fiat-plate  model  in  the  CAL  48"  Hypersonic 
Shock  Tunnel  to  completely  span  the  transition  between  the  classical  thin  bound¬ 
ary  layer  regime  and  the  strong  interaction  regime.  Both  of  these  facilities 

3  7  38 

have  been  reported  in  available  literature  *  and  only  those  aspects  of  the 
facilities  that  have  direct  bearing  on  the  present  results  will  be  desc ribed  here. 

A.  The  Six-Foot  Hypersonic  Shock  Tunnel 

The  CAL  Six-Foot  Shock  Tunnel  utilizes  the  tailored-interface  principle  of 
operation.  Its  important  features,  with  regard  to  the  present  research,  are 
the  wide  range  of  test  conditions  that  can  be  obtained,  the  mode  of  operation, 
and  the  uniformity  of  the  flow  in  the  test  section. 

The  Six-Foot  Shock  Tunnel  is  designed  to  produce  reservoir  pressures  of 
25-2000  atmospheres  and  reservoir  temperatures  of  4000  -8000°K  in  uncontam¬ 
inated  air.  The  tunnel  is  currently  operated  at  reservoir  pressures  up  to  1300 
atm.  and  at  reservoir  temperatures  of  4000  -  7000°K.  The  present  investigations 
were  made  at  a  nominal  reservoir  temperature  of  about  4000"K,  and  over  a 
pressure  range  of  25  -  450  atmospheres.  The  mode  of  operation  selected  for  the 
research  was  to  use  a  heated  driver  (675°K)  consisting  of  a  mixture  of  two  parts 
of  hydrogen  and  one  part  of  helium.  With  this  mixture,  effective  interface  tai¬ 
loring  was  obtained  at  a  shock  Mach  number  of  about  6.  25,  producing  a  reservoir 
temperature  of  about  4000,K.  This  particular  reservoir  temperature  was 
selected  tu  provide  a  teBt  section  static  temperature  high  enough  to  prevent 
liquefaction  and  to  insure  that  the  thermochemical  nonequilibrium  effects  in  the 
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nozzle  expansion  would  be  small.  This  latter  consideration  is  discussed  more 
fully  in  the  next  section. 

B.  Tunnel  Calibrations 

An  important  consideration  in  low-density  experiments  is  that  of  flow  uni¬ 
formity  in  the  test  section.  The  flows  were  calibrated  in  the  present  investi¬ 
gation  at  each  test  condition.  The  calibrations  consisted  of  measuring  the 

distribution  of  pitot  pressures  across  the  test  section,  using  piezoelectric  pres- 
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sure  transducers  developed  at  CAL,  and  measuring  the  flow  angularity  on 
the  tunnel  centerline  with  a  pair  of  flat-plate  models. 

The  pitot  pressure  calibrations,  Fig.  1,  demonstrate  that  at  low  densi¬ 
ties,  the  wall  boundary  layer  tends  to  close.  The  lowest  test  section  density 
(Re  s*  350  ^ /in.  )  produced  a  useable  inviscid  core  of  about  6  inches.  At  the 
higher  reservoir  pressures,  the  flow  is  acceptably  uniform  over  at  least  the 
center  3  feet  of  the  test  section,  and  there  is  no  evidence  of  the  wall  boundary 
layer. 

The  flow  angularity  calibration  was  based  on  the  experimental  observation 
that  at  low  densities,  a  region  of  nearly  constant  heat  transfer  exists  over  a 
considerable  portion  of  a  sharp  flat  plate  near  the  leading  edge,  A  pair  of  small 
scale  models  (Fig.  2)  were  built  to  obtain  detailed  heat  transfer  information  in 
this  region.  These  were  mounted  with  the  instrumented  surfaces  parallel  and 
facing  each  other,  and  tested  in  a  null -type  experiment.  The  angle  of  attack  of 
the  array  was  varied  over  a  small  range  to  determine  the  flow  angularity;  that 
is,  the  angle  at  which  both  models  indicated  the  same  heat  transfer  rates.  The 
advantage  in  this  technique  is  that  in  the  approach  to  free-molecule  flow,  the 
heat  transfer  on  a  flat  plate  is  remarkably  sensitive  to  angle  of  attack.  Typically, 
a  one  degree  change  in  angle  of  attack  produces  a  70%  change  in  heat  transfer 
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rate.  Consequently,  flow  angularity  is  measured  with  a  precision  of  about 

+  0.  1*  . 

The  heat  transfer  rate  used  in  these  calibrations  was  an  average  value  in 
the  region  where  the  heat  transfer  rate  was  essentially  invariant  with  the  chord- 
wise  position.  Typical  results  obtained  at  the  lowest  density  condition  used  in 
this  research  are  presented  in  Fig.  3.  The  data  in  the  upper  half  of  the  figure 
are  plotted  as  a  function  of  geometric  angle  of  attack,  and  indicate  a  flow 
angularity  of  about  0.  1°.  The  self-consistency  cf  these  data  is  demonstrated 
in  the  lower  half  of  the  figure  where  the  indicated  angularity  has  been  applied 
to  the  data  from  both  models,  and  plotted  as  a  function  of  aerodynamic  angle  of 
attack.  The  fact  that  the  data  from  both  models  produce  a  smooth  curve  demon¬ 
strates  that  the  inferred  flow  angularity  is  valid.  The  measured  angularity  at 
the  other  test  conditions  was  0.  1  -  0.  15°  ;  it  is  concluded  that  the  effects  of  flow 
angularity  are  comparable  with  or  smaller  than  the  experimental  scatter. 

C.  Models  and  Instrumentation 

Two  flat-plate  models  were  used  in  the  experimental  research.  The  small- 
scale  model,  shown  in  Fig.  2,  was  designed  to  investigate  the  leading-edge 
region,  and  consisted  of  fourteen  thin-film  (platinum)  thermometers  mounted  on 
a  Pyrex-brand  glass  plate.  This  glass  plate  was  bonded  to  a  steel  plate. 
Typically,  these  thermometers  were  ,  02"  x  .  25"  and  near  the  leading  edge  were 
spaced  about  .  02"  apart.  All  fourteen  gages  were  located  within  about  1"  from 
the  leading  edge,  and  under  typical  test  conditions,  five  gages  were  within  one 
ambient  mean-free-path  from  the  leading  edge.  The  first  gage  was  at  the  lead¬ 
ing  edge.  The  gage  lead  wires  were  gold  films  about  2  microns  thick,  leading 
to  conventional  lead  wires  at  the  model  extremities.  The  entire  model  was 
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coated  with  titanium  dioxide  about  0.  1  microns  thick  to  prevent  electrical  short¬ 
ing  in  ionized  flows.  The  operating  principles  of  the  thin-film  gage  as  well  as 

40  41  42 

the  coated  gage  have  been  described  in  the  literature,  ’  ’  and  will  not  be 

reviewed  here. 

The  leading  edge  was  flat-faced  with  a  thickness  estimated  to  be  less  than 
.  005".  The  wedge  angle  of  the  lower  surface  was  30°.  During  the  experiments 
at  zero  angle  of  attack,  it  was  found  that  the  first  gage  indicated  a  heat  transfer 
rate  about  twice  that  at  the  next  station.  It  has  been  determined  that  this 
behavior  was  due  to  heat  conduction  from  the  blunt  leading  edge  and  the  lower 
surface.  Under  f ree-molecule  conditions,  the  heat  transfer  to  the  blunt  leading 
edge  and  to  the  wedge  lower  surface  is  respectively  sixty  times  and  thirty  times 
greater  than  that  to  the  upper  surface.  Heat  conduction  calculations  showed 
this  was  sufficient  to  cause  a  100%  change  in  indicated  heat  transfer  to  the  in¬ 
strumented  surface.  Experiments  were  made  with  the  instrumented  surface 
shielded  from  the  flow  to  measure  the  heating  from  both  surfaces.  These  veri¬ 
fied  the  above  calculation.  .This  is  the  first  instance  known  to  the  authors  where 
the  assumption  of  one -dimensional  heat  transfer  was  not  valid.  Careful  con¬ 
sideration  of  these  effects  are  necessary  when  investigating  low-density 
leading-edge  phenomena.  Those  data  that  were  affected  by  this  additional  heat¬ 
ing  are  not  reported. 

The  large-scale  flat  plate  consisted  of  a  steel  plate  with  detachable  lead¬ 
ing  and  trailing  edges.  The  model  inst  rumentation  shown  in  Fig.  4  included 
platinum  thin-film  resistance  thermometers  and  piezoelectric  pressure  trans¬ 
ducers.  The  primary  thermometers  and  pressure  transducers  were  equi¬ 
distant  from  the  midspan  position.  Thermometers  at  spanwise  locations  were 
included  to  provide  a  continuous  check  on  finite  span  effects.  Similarly,  three 
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chordwise  extensions  were  used  to  obtain  a  quantitative  measure  of  effects 
associated  with  the  finite  chord  length.  The  thermometers  were  on  glass  but¬ 
tons,  mounted  in  removable  sleeves  to  facilitate  instrumentation  changes.  The 
detachable  leading  edge  was  instrumented  with  thermometers  up  to  0.  3  inches 
from  the  leading  edge,  and  with  a  pressure  transducer  0.  61  inches  from  the 

_4 

leading  edge.  The  leading  edge  was  honed  to  a  thickness  of  the  order  of  5  x  10 
inches. 

The  piezoelectric  transducers  are  those  described  in  Ref.  39.  Briefly, 

they  are  diaphragm-type  transducers  driving  modified  lead-zirconium-titanate 

crystals.  The  nominal  output  is  2  volts/psi,  and  the  transducers  have  been 

used  in  the  present  research  to  measure  pressures  as  low  as  2  x  10"^  psi.  An 

important  consideration  in  low-density  experiments  is  the  lag  time  associated 

with  pressure  orifice  size  and  tube  internal  volume.  This  aspect  was  examined 

43 

within  the  framework  of  both  continuum  and  free-molecule  concepts.  It  was 

-3  2 

concluded  that  with  an  orifice  area  of  1  5  x  10  in.  ,  an  orifice  length  of 

.  2  -  3  3 

6x10  in.  ,  and  an  internal  volume  of  11  x  10  in.  ,  the  lag  time  should  be 

about  200  microseconds.  In  the  subsequent  experiments  it  was  impossible  to 

quantitatively  assess  the  lag  time  since  it  was  smaller  than  the  tunnel  starting 

time  (about  1  millisecond).  However,  the  pressure  data  were  found  to  achieve 

a  steady  value  within  200  -  300  microseconds  after  the  thermometers  showed  the 

tunnel  had  started. 

In  the  interests  of  aerodynamic  cleanliness,  the  large-scale  model  was 
supported  by  8  wires  1/16"  in  diameter  and  attached  at  the  model  extremities. 
The  instrumentation  lead  wires  were  let  out  of  the  model  under  surface  and 
hidden  in  the  shadow  of  the  lower  surface  support  wires.  The  support  wires 
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were  connected  to  the  tunnel  walls  through  springs  to  minimize  accelerations 
that  could  excite  the  pressure  transducers. 

All  data  were  recorded  on  standard  oscilloscopes  and  high-speed  oscillo¬ 
graphs.  Special  instrumentation  was  constructed  at  CAL  to  record  the  low  heat 
transfer  levels  measured.  In  order  to  eliminate  the  numerical  data  reduction 
procedures  to  relate  the  measured  surface  temperatures  to  heat  transfer  rates, 

analogue  circuits  have  been  devised  to  accomplish  the  numerical  integration 

44 

during  the  experiment.  However,  these  are  effectively  filters  and  thereby 
attenuate  the  electrical  signal.  Recent  modifications  were  made  in  the  circuit 
networks  to  include  low-noise  amplifying  stages  which  decreased  the  noise  level 
by  about  two  orders  of  magnitude. ^  During  related  research,^1®  heat  transfer 
levels  as  low  as  2  x  10"^  Btu/ft^  sec.  (corresponding  to  10"^  *K  temperature 
change  in  1  millisecond)  were  measured  with  a  signal -to -noise  ratio  of  about  15, 
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IV.  AMBIENT  TEST  SECTION  CONDITIONS 


The  experiments  in  the  48-inch  Hypersonic  Shock  Tunnel  were  made  at 
reservoir  temperatures  of  about  3000°K  and  at  reservoir  pressures  of  50-400 
atm.  These  conditions  are  sufficiently  moderate  that  in  the  nozzle  expansion, 
the  air  can  be  regarded  as  in  the rmo-chemical  equilibrium.  Consequently,  the 
ambient  test  section  conditions  were  determined  following  the  equilibrium  cal¬ 
culation  described  below. 

As  noted  earlier,  the  low  density  experiments  in  the  Six-Foot  Hypersonic 
Shock  Tunnel  were  made  at  a  reservoir  temperature  of  about  4000°K  and  with 
reservoir  pressures  of  25-450  atm.  These  conditions  were  selected  to  prevent 
liquefaction  in  the  test  section  and  to  minimize  the  thermochemical  nonequilibrium 
effects  in  the  nozzle  expansion.  The  alterations  due  to  nonequilibrium  were  taken 
into  account  by  calculating  the  test  section  conditions  assuming  an  equilibrium 
nozzle  expansion,  and  then  introducing  the  changes  stemming  from  nonequilibrium. 

The  reservoir  conditions  for  near-tailored-interface  operation  are  deter¬ 
mined,  using  equilibrium  air  calculations,  by  measuring  the  speed  of  the  incident 
shockwave  and  the  reservoir  pressure.  The  departures  from  ideal  tailored- 
interface  conditions  are  taken  into  account  by  assuming  isentropic  wave  processes 
caused  the  observed  departures.  The  equilibrium  test  section  conditions  are 
determined,  assuming  an  isentropic  expansion,  from  the  measured  pitot  pressure 
and  the  hypersonic  approximations 


where  £  is  the  density  ratio  across  a  normal  shock  wave  and  to  is  the 
pitot  pressure.  The  results  of  this  procedure  have  been  checked  with  exact 
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numerical  results  and  have  been  found  to  be  accurate  to  +  0.  4%  over  the  range 
of  test  conditions  covered  here. 

The  effects  of  thermochemical  nonequilibrium  in  the  nozzle  expansion  were 

45 

calculated  using  the  computer  program  developed  by  Eschenroeder,  et  al. 

Results  obtained  with  this  machine  program  were  recently  compared  with  limited 

46 

experimental  data  by  Hall,  et  al  and  found  to  be  in  good  agreement.  Equilibrium 
and  nonequilibrium  solutions  were  obtained  for  reservoir  conditions  approximating 
those  obtained  in  the  experiments;  To  =  4000-K,  Fo  =  25,  100,  300  atm. 

These  solutions  were  run  only  to  an  axial  station  in  the  nozzle  expansion  where 
the  flow  was  completely  frozen  and  could  be  regarded  as  ideal.  The  solutions  were 
then  extended  to  the  test  section  conditions  by  assuming  an  isentropic  expansion 
and  a  specific  heat  ratio  ^  =1.4.  The  important  item  to  be  noted  in  this  pro¬ 

cedure  is  that  both  the  equilibrium  and  nonequilibrium  solutions  must  be  extended 
to  the  measured  pitot  pressure. 

A  discussion  of  the  aerodynamic  parameters  pertinent  to  the  present  research 
is  given  in  Sections  II  and  V.  A  tabulation  of  the  nonequilibrium  corrections  to  these 
parameters  is  given  in  Table  I  for  each  of  the  three  nominal  reservoir  conditions 
used  in  this  research.  The  subscripts,  ft  and  6  ,  refer  to  nonequilibrium 

and  equilibrium  values,  respectively.  The  important  item  to  be  noted  in  Table  I 
is  the  magnitude  of  the  nonequilibrium  corrections  to  certain  of  the  established  aero¬ 
dynamic  parameters.  For  example,  consider  the  strong  interaction  parameters, 

M  Ch  ,  — and  ^  ,  at  a  reservoir  pressure  of  25  atm.  The  non- 

leo 

equilibrium  corrections  are  respectively  53%,  26%,  and  14%.  This  illustrates 
that  careful  attention  must  be  given  to  nozzle  nonequilibrium  effects  when  inter¬ 
preting  conventional  aerodynamic  test  information  obtained  even  at  the  modest  tem¬ 
perature  of  4000  “K. 
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TABLE  I 


Nozzle  Nonequilibrium  Effects  on  the  Governing 
Aerodynamic  Parameters 

”To  =  4000°K 
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V.  EXPERIMENTAL  RESULTS 

The  experimental  program  consisted  of  testing  the  large-scale  and  small- 
scale  flat  plates  at  zero-angle  of  attack  and  at  large  angles  of  attack.  Surface 
pressure  and  heat  transfer  rates  were  measured  in  all  experiments  with  the 
large-scale  flat  plate,  but  only  heat  transfer  data  will  be  reported  for  the  wedge- 
flow  experiments.  The  pressure  data  from  those  experiments  generally  confirmed 
the  theoretical  conclusion  that  the  surface  pressure  would  be  the  ideal  wedge  pres¬ 
sure.  The  wedge  flow  experiments  will  be  discussed  first  to  demonstrate  the 
dominant  effects  due  to  transport  processes  at  the  shock  wave. 

A.  Wedge  Flow  Results 

As  indicated  previously,  one  important  aspect  of  low  density  wedge  flows  is 

that  they  offer  a  means  for  isolating  a  single  fluid  dynamic  mechanism  and  studying 

17 

it  in  detail.  In  the  theoretical  study  it  was  shown  th;  „  the  effects  of  all  other  mech¬ 
anisms  were  numerically  small  in  comparison  with  that  due  to  the  transport  processes. 
The  purpose  of  the  experiments  is  to  verify  the  theoretical  results,  duplicating  as 
closely  as  possible  the  theoretical  assumptions. 

One  important  experimental  requirement  in  checking  the  theory  is  on  the  wedge 
angle.  The  competing  low  density  effects  can  be  suppressed  and  certain  of  the  the- 


oietical  assumptions  are  valid  only  if  the  wedge  angle  is  sufficiently  large.  However, 
the  theory  assumes  that  the  shock  wave  is  attached,  which  requires  that  the  wedge 
angle  be  less  than  about  42°.  For  these  reasons,  the  verifying  experiments  were 
made  with  wedge  angles  of  20°,  30°,  and  40°.  These  results  are  presented  and  com¬ 
pared  with  the  theory  in  Fig.  5.  It  can  be  seen  that  the  experimental  results  are 
generally  in  good  agreement  with  the  theory  with  the  experimental  data  falling  some¬ 


what  above  the  theory.  In  the  high  density  limit  (large  values  of 
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the  theory  approaches  the  Blaeius  solution  for  |r  =  1.  The  average  of  the  exp¬ 
erimental  data  falls  some  10%-15%  above  the  Blasius  limit,  suggesting  that  the 
difference  between  theory  and  experiment  can  be  attributed  to  the  assumption  of 
Ft  =  1.  As  a  check  on  this  hypothesis,  a  semi-empirical  formula^  which  takes 
account  of  Prandtl  numbers  different  from  unity  was  compared  with  the  results  in 
the  high  density  limit.  This  approximate  calculation  was  about  7%  above  the  theory, 
lending  support  to  the  hypothesis  that  the  observed  differences  between  theory  and 
experiment  are  due  to  differences  in  Prandtl  number. 

An  important  consideration  in  wedge  flow  experiments  is  the  effect  of  a 
finite  wedge  length.  The  theory  is  based  upon  a  semi-infinite  model,  while  exper¬ 
iments  must  use  finite  model  lengths.  The  experiment  then  has  a  9trong  expansion 
fan  situated  at  the  model  trailing  edge  which  can  influence  the  upstream  portion  of 
the  model  through  the  subsonic  portions  of  the  viscous  layer.  This  effect  can  be 
evaluated  by  comparing  the  data  obtained  with  the  large  scale  and  small  scale  models 
at  comparable  test  conditions.  Some  of  the  instrumentation  on  both  models  were 
at  equivalent  chordwise  stations.  The  effects  of  finite  chord  length  should  be  neg¬ 
ligibly  small  for  the  leading  edge  instrumentation  on  the  large  model.  If  there  is 
a  trailing  edge  effect,  it  should  be  apparent  in  the  data  from  the  small  scale  model. 
The  comparison  in  Fig.  6  shows  the  data  from  the  two  models  agree  well  in  the 
regions  where  there  was  instrumentation  duplication  and  that  there  is  no  system¬ 
atic  difference  between  the  data  from  the  two  models. 

The  range  of  wedge-angles  over  whirh  the  theory  is  applicable  is  important, 
The  data  reported  here  show  the  theory  is  valid  for  wedge  angles  of  20o-40o;  exper¬ 
iments  are  now  in  progress  to  better  define  that  range.  Preliminary  results  show 
that  the  theory  is  not  applicable  for  wedge  angles  of  10",  and  is  not  applicable  near 
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the  leading  edge  for  wedge  angles  that  cause  shock  wave  detachment.  The 
experiments  indicate  the  theory  does  apply  over  certain  portions  of  the  wedge 
for  angles  as  large  as  75°. 

Two  conclusions  from  the  wedge  flow  results  have  bearing  on  the  zero- 
angle  of  attack  experiments.  First,  the  theory  is  adequate  over  the  entire  density 
range  (from  the  thin  boundary  layer  regime  to  the  free-molecule  limit),  although 
small  inaccuracies  (10%-15%)  apparently  stem  from  the  assumption  of  Rr  =  1. 
Second,  since  the  dominant  fluid  mechanism  causing  the  low  density  departures 
from  the  thin  boundary  layer  regime  is  the  transport  processes,  these  results 
provide  a  basis  for  determining  the  conditions  when  transport  effects  first  become 
evident.  Comparing  the  faired  curve  (averaging  the  high  density  data)  with  the  low 
density  data,  it  can  be  seen  that  the  low  density  departures  begin  at 

-  1  ...  i&r  .3  to  .4.  In  terms  of  the  parameter  to  be  used  in  examining 

cosoc  i  MzCx 

the  zero  angle  of  attack  data,  this  expression  conservatively  corresponds  to 

M/S;*  i  • 

B.  Zero-Angle  of  Attack 

The  flat  plate  experiments  at  zero  angle  of  attack  covered  a  range  of  Mach 
numbers  and  Reynolds  numbers  to  include  the  transition  from  the  weak  to  the  strong 
interaction  regime  through  the  approach  to  free-molecule  flow.  Heat  transfer 


measurements  will  be  reported  for  Knudsen  numbers  as  large  as  — ~  =  5, 

and  pressure  measurements  will  be  reported  for  Knudsen  numbers  as  large  as 
1/5.  Since  the  greatest  body  of  theoretical  work  is  derived  from  the  continuum  view¬ 
point,  the  data  will  be  discussed  with  the  aim  of  improving  the  understanding  in  this 
area  and  explaining  the  departures  from  established  continuum  theory.  Separate 
comparisons  will  be  made  with  the  applicable  leading  edge  solutions. 
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Transition  Between  Weak  and  Strong  Interactions.  The  heat  transfer 
and  pressure  data  obtained  at  zero  angle  of  attack  have  been  cast  into  the  viscous 
interaction  parameters  and  are  presented  in  Figs.  7  and  8.  The  heat  transfer 
data  are  compared  with  the  theory  for  zero  pressure  gradient,  the  weak  inter¬ 
action  theory  derived  from  the  first-order  result  for  the  pressure  (Eq.  3),  and 

4 

with  strong  interaction  theory  from  Ref.  9  (Eq.  5).  Lees'  theory  is  nearly 
coincident  with  that  shown,  and  is  not  included  here.  The  pressure  data  are 
compared  with  the  first-order  weak  and  the  strong  interaction  theories,  Eqs. 

(1)  and  (4),  as  well  as  with  that  of  Ref.  4.  It  can  be  seen  that  there  are  notable 
departures  from  strong  interaction  theory  at  small  and  large  values  of  ^ 

The  departures  exhibited  by  the  heat  transfer  data  at  low  values  of  ^ 

define  the  transition  to  the  weak  interaction  regime.  This  behavior  was  critically 

9  13 

examined  in  previous  experiments.  ’  This  transition  in  heat  transfer  is  accurately 
predicted  within  weak  interaction  theory  if  the  higher  order  terms  are  retained  in 
the  heat  transfer  formula,  Eq.  (3).  Figure  7  indicates  that  the  transition  between 
the  two  regimes  is  rather  narrow,  and  for  practical  purposes,  it  is  defined  by 
5  (  ^  (,  20.  The  correlation  of  the  pressure  data  with  weak  and  strong  inter¬ 

action  theories,  Fig.  8,  shows  larger  departures.  The  data  at  the  smallest  values 
of  ^  agree  well  with  both  weak  and  strong  interaction  theory,  indicating  the  tran¬ 
sition  regime. 

Strong  Shock  Wave  Approximations.  it  was  noted  earlier  that  the  strong 
shock  wave  approximations  evidently  could  impose  severe  limitations  on  strong 
interaction  theory,  especially  as  they  influence  the  density  ratio  The  inferred 
limits,  given  by  Eqs.  (8)  and  (9),  suggest  the  theory  should  break  down  for 


23.  The  agreement  between  theory  and  experiment  for 
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10-20,  especially  away  from  the  leading  edge  (decreasing  ^  ), 

suggest  that  the  strong  shock  wave  approximations  do  not  impose  any  effective 
limit  on  theory.  Indeed,  the  agreement  between  theory  and  experiment  is  as 
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The  data  in  Figs.  7  and  8  show  large  departures  from  strong  interaction 
theory  at  large  values  of  ^  .  This  behavior  has  been  experimentally  ob¬ 
served  by  others.  ^  ^  In  that  work,  it  was  hypothesized  that  these  departures 

stemmed  from  the  effects  of  velocity  slip  and  temperature  jump  at  the  surface. 

This  conclusion  was  reached  on  the  basis  of  correlations  similar  to  those  in  Figs. 

7  and  8  and  on  the  basis  of  Schlieren  observations.  In  view  of  the  competing  low 
density  effects  described  in  Section  II  that  exist  in  the  strong  interaction  regime, 
and  the  importance  of  identifying  the  governing  mechanisms,  a  more  critical 
assessment  of  this  hypothesis  is  warranted. 

Low  Density  Effects  on  Viscous  Interactions.  Before  proceeding  with  the 
discussion  of  the  higher-order  fluid  dynamic  mechanisms,  it  must  be  noted  that 
some  of  the  low  density  departures  observed  in  Figs.  7  and  8  could  stem  from 
finite  trailing  edge  effects.  To  check  this  possibility,  the  large  scale  flat-plate 
model  was  tested  at  one  of  the  lower  density  conditions  with  three  different  trailing 
edge  configurations;  the  test  trailing  edge  extension  (5"  long),  no  trailing  edge 
extension,  and  a  10"  trailing  edge  extension.  Typical  results  of  this  investigation 
are  shown  in  Fig.  9.  The  form  for  presenting  the  data  here  will  be  developed  sub¬ 
sequently.  The  important  item  is  the  fact  that  the  data  with  all  three  trailing  edge 
configurations  agree  within  the  scatter  of  the  data,  demonstrating  that  no  discernable 
effect  due  to  a  finite  trailing  edge  was  present. 

It  was  stated  earlier  that  four  effects  neglected  in  the  strong  interaction 
regime  are  the  vorticity  interaction,  the  shockwave  heating,  the  effects  of  surface 


35 


slip  and  the  effects  of  transport  processes  at  the  shockwave.  Moreover,  these 

M/S’ 


four  effects  are  all  governed  by  the  same  parameter, 


Re, 


and  because 


of  the  perturbation  analysis  used  in  the  derivations,  the  effects  all  enter  in  the 
same  manner.  - - ,  ^ 

-|r  ~  %  [  1  +  const: 

M3C«  ~  X^[l  +  CONSX(M/^) 


The  importance  of  these  four  effects  in  the  strong  interaction  regime  can  then  be 

-1-EL  v.  M\/jP^  and  M^Ch 

%  p-  v 


examined  in  correlations  of 

T5T 


m/ 


Re, 


These  correlations  are  shown  in  Figs.  10  and  11  and 


are  compared  with  strong  and  weak  interaction  theories.  The  equations  for  each 


of  these  theories  in  terms  of 
Weak  Interaction: 


follow  from  Eqs.  (1),  (3),  (4),  and  (5). 
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+- 
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where 


Strong  Inte 


ft  =  A(Prj^-2B(Pr)] 
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Strong  interaction  theory  appears  in  Figs.  10  and  11  as  a  horizontal  line 
which  varies  only  with  the  wall  temperature  ratio.  The  departures  from  the 
strong  interaction  regime  are  then  evident  as  an  increase  and  a  decrease  in 


at  the  weak  interaction  limit  and  the  low  density 


X  P  and  M3Ch 
t  %  i- 

limit,  respectively.  It  should  be  noted  in  Eqs.  (34)  and  (35)  that  this  is  not  a 


generalized  correlation  in  that  there  is  a  Mach  number  dependence  in  the  weak 
interaction  regime. 

Figure  10  shows  that  the  heat  transfer  data  correlate  well  in  both  the  weak 
interaction  and  the  strong  interaction  regimes,  and  that  this  correlation  is  quite 
satisfactory  for  defining  a  low-density  limit  to  the  strong  interaction  regime; 


that  is,  M 


^  1/5.  The  scatter  in  the  weak  interaction  and  strong  inter¬ 


action  regimes  is  about  +  8%,  typical  of  experimental  accuracy  in  shock  tunnels. 
For  values  of  Mv/gs-  >  1/5  ,  the  data  in  Fig.  10  are  seen  to  scatter  and  to 

Y  Ken 

separate  into  two  fairly  distinct  branches,  corresponding  to  the  higher  and  lower 
density  test  conditions.  Examination  of  Fig.  8  shows  that  for  the  two  runs  at  the 
lowest  densities,  the  pressure  was  essentially  constant  over  the  enLire  model, 
analogous  to  that  in  free-molecule  flow.  In  contrast,  the  pressure  variation  for 
the  two  higher  density  runs  was  about  half  of  that  predicted  by  strong  interaction 
theory.  It  was  noted  earlier  the  parameter  governing  the  vorticity  interaction, 

k  A  fC% 

the  shock  wave  heating,  and  the  slip  effects  was  |'v/)  y' only  if  the  pressure 
distribution  was  that  given  by  strong  interaction  theory.  The  lack  of  correlation 
for  >  ^1/5  can  be  interpreted  as  a  consequence  of  the  fact  that  the 

pressure  is  no  longer  adequately  described  by  strong  interaction  theory. 

The  pressure  data,  Fig.  11,  qualitatively  show  the  same  behavior  as  the 
heat  transfer  data,  except  that  the  departures  are  more  evident.  There  is 
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virtually  no  strong  interaction  regime  under  the  test  conditions  for  the  present 
experiments.  Indeed,  these  suggest  the  limit  on  the  strong  interaction  regime 

As  noted  earlier,  experimentally  it  is  not  feasible  to  distinguish  between 
the  four  effects  encountered  in  the  strong  interaction  regime  because  all  are 
essentially  of  the  same  order.  It  can  be  concluded,  however,  that  the  net  result 
of  all  four  effects,  in  the  case  of  the  cold  wall,  is  to  cause  a  decrease  in  pres- 
sure  and  heat  transfer  for  values  of  Pv|y|=j —  ^  1/10  -  1/5  . 

Generalization  of  Governing  Parameters.  It  is  evident  that  the  separate 
effects  of  the  vorticity  interaction,  shock  wave  heating,  transport  processes  at 
the  shock,  and  the  slip  boundary  conditions  cannot  be  examined  within  the  frame¬ 
work  of  strong  interaction  theory.  But  it  is  important  that  the  mechanism  governing 
the  observed  departures  from  strong  interaction  theory  be  better  identified  to  guide 
further  theoretical  developments.  Considerable  insight  into  the  dominant  mech¬ 
anism  can  be  obtained  by  generalizing  the  parameter  governing  each  of  the  higher  - 
order  effects.  It  was  previously  observed  that  the  data  correlation  in  N^l  y 
might  cease  to  be  valid  because  the  strong  interaction  description  of  the  pressure 
distribution  becomes  inaccurate.  The  generalization  of  the  governing  parameters 


then  consists  of  retaining  the  pressure  terms,  as  opposed  to  assuming. 


% 


Neglecting  the  details  of  the  development,  it  can  be  shown  that  each  of  the  higher- 
order  effects  is  governed  by  the  following  parameters. 


Vorticity  Interaction 


f 


|  _  _ 

1^.  F ll 


IE]1 


Shock  Wave  Heating 


(l&)  ~  ,f 


4-  /  57  \213T 
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(38) 


(39) 
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Surface  Slip  and  Shock  Transport  Processes 


u) 


M 


a 


r-s^j 


(40) 


The  dominant  mechanisms  can  now  be  identified  by  inspection.  Consider 
the  heat  transfer  data  in  Fig.  10  and  the  effect  of  substituting  each  of  the  para¬ 


meters,  Eqs.  (30),  (39  ),  and  (40  ),  for 


[7^p 

Mvi?S7  .  The  pressure  data. 
Fig.  12,  shows  that  both  the  vorticity  interaction  parameter  and  the  shock  wave 
heating  parameters ,  when  substituted  for  |^J  > 


would  tend  to  aggravate 


the  correlation  of  the  heat  transfer  data.  The  pressure  level  under  the  lowest 
density  conditions  is  smaller  than  that  at  higher  densities.  Consequently,  when 
substituting  Eqs.  (38  )  and  (  39)  for  ,  the  two  low  density  branches 

would  tend  to  separate.  In  contrast,  Eq.  (40)  shows  that  parameter  varies  inver¬ 
sely  with  the  square-root  of  the  pressure  ratio,  which  would  tend  to  bring  the 
low-density  data  into  agreement.  Hence,  an  improved  data  correlation  should 

i3r\.  ip 

and  __  1  as  a  function  of  the 

parameter  in  Eq.  (40). 

One  can  employ  a  semi-empirical  procedure  to  show  that  a  consistent  data 


be  obtained  in  terms  of  M 

*♦ 


correlation  can  be  obtained  in  terms  of  the  parameter  governing  surface  slip  and 
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the  transport  effects.  This  was  done  by  Vidal  and  Wittliff,  and  consists  of 
noting  in  Fig.  12  that  the  pressure  distribution  for  each  experiment  is  well  rep¬ 


resented  by 


£.  „_L  . 


where  the  exponent,  K)  ,  varies  with  ambient 


density.  The  parameter  in  Eq.  (40  )  can  then  be  integrated  to  yield 

-  I^C«(i-n)l1: 

J 

i  can  then  be 


(41) 


The  heat  transfer  and  pressure  data 


correlated  consistently  in  plots 
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of 


M3Ch 


and 


1  £_  v 

7-  Be  V*  ^ 


Such  a  correlation  is  of  practical 


%T  r  loo 

use  in  that  it  graphically  demonstrates  that  the  slip  effects  and/or  transport 

effects  dominate  over  the  vorticity  interaction  and  shock  wave  heating  effects 
in  this  low-density  regime.  The  question  still  remains  as  to  the  relative  impor¬ 
tance  of  surface  slip  and  the  transport  processes  at  the  shock  wave. 

Surface  Slip  and  Shock  Transport  Processes.  The  relative  magnitudes 
of  surface  slip  and  the  effects  of  shock  transport  processes  in  the  present  experi¬ 
ments  can  be  assessed  by  comparing  the  boundary  conditions  for  each,  Eqs.  (10), 
(11),  (  14  ),  and  (  17).  The  ratio  of  the  velocity  decrement  at  the  shockwave 

to  the  slip  velocity  at  the  surface  and  the  ratio  of  the  enthalpy  decrement  at  the 


shock  wave  to  the  enthalpy  jump  at  the  surface  are  approximately 

Us  *  .  x.  (M 

I  4  ,  and  that  the  measured  sur- 
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Noting  that  for  low  temperatures, 


(42) 


face  pressure  can  be  used  in  the  oblique  shock  relations  to  calculate  the  temperature 
ratio  across  the  shock  wave, 

Ik  *  ZYCX-fl  +  _4.y  „  -  zOf-i) 

Too  tthO1  U-h;  u+il5 


tvpf*  *  mu 

T  2K  P* 


it  is  possible  to  estimate  the  relative  magnitudes  of  the  slip  effects  and  the  shock 
transport  effects.  The  ambient  stream  temperature  in  these  experiments  was 


about  60°K,  and  the  wall  temperature  was  about  300°K.  The  ratio  of  the  mean- 

\s 

\w 


free-path  at  the  shock  wave  to  that  at  the  surface  was 


3/4  to  1/8, 


depending  on  the  ambient  density  level.  The  larger  value  of  this  ratio  occurred 
during  the  higher  density  experiments  where  the  departures  from  strong  interaction 
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theory  were  small.  If  this  is  to  be  ascribed  to  the  shock  layer  becoming  fully 
viscous,  this  experiment  would  then  correspond  to  the  initial  oscillation  between 
the  shock  wave  and  the  viscous  layer.  The  velocity  gradient  at  the  shock  wave 
would  then  be  small  in  comparison  with  that  at  the  surface,  and 


At  the  lower  densities, 


1/8,  and  if  the  velocity  gradient  at  the 


7\w 

shock  is  taken  to  be  comparable  to  that  at  the  wall,  the  effects  of  transport  pro¬ 
cesses  at  the  shock  wave  would  still  be  an  order  of  magnitude  smaller  than  the 
effects  of  surface  slip.  Consequently,  it  is  concluded  that  the  effects  of  trans¬ 
port  processes  at  the  shockwave  were  always  small  in  comparison  with  those 
due  to  surface  slip  in  the  present  experiments, 

This  conclusion  is  verified  by  the  experimental  data.  In  the  discussion  of 
wedge  flows  it  was  noted  that  shock  wave  slip  was  the  dominant  mechanism  there 


and  that  the  effects  of  shock  wave  slip  were  first  apparent 


fo,  M'/U 


»  1.0. 


Examination  of  Figs,  10  and  11  shows  that  the  low-density  departures  become 
important  at  much  higher  densities  (M\ZS!a  1/  10  -  1/^than  in  the  wedge - 
flow  case.  Moreover,  the  temperature  at  the  shock  wave  in  the  wedge-flow  case 
is  about  an  order  of  magnitude  greater  than  in  the  zero-angle  cf  attack  case.  For 
these  reasons,  it  is  concluded  that  the  dominant  mechanism  causing  the  observed 
low-density  departures  is  slip  at  the  surface. 

Effects  of  the  Pressure  Jump.  The  final  argument  to  demonstrate  that  the 
observed  departures  from  strong  interaction  theory  stem  from  surface  slip  is  to 
compare  these  departures  with  the  pressure  jump.  As  noted  earlier,  the  theory  of 
Ref.  30  is  based  on  a  perturbation  method  applied  to  strong  interaction  theory  to 
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examine  the  first-order  slip  effects.  The  theory  is  applicable  only  for  the  case 
of  an  adiabatic  wall,  due  to  the  approximations  invoked  to  effect  a  solution.  In 
addition,  this  theory  treats  the  classical  velocity  and  temperature  jump  effects 
but  neglects  the  pressure  jump.  The  pressure  data  are  compared  with  the 
theory^  in  Fig.  11.  The  pressure  jump  as  given  earlier  (Eq,  24}  and  as  given 
in  Ref.  29  has  been  applied  to  the  results  of  the  slip  analysis  as  follows, 


fit  = 


(43) 


j_  FLi  r&i 
[%  fiJ  LRj 

where  the  first  term  on  the  right  is  the  result  given  in  Ref.  30  and  reproduced 
here  as  Eq.  29  ,  and  the  second  term  is  the  pressure  jump.  It  can  be  seen  that 
the  changes  in  the  induced  pressure  due  to  alterations  in  the  boundary  layer  dis¬ 
placement  thickness  are  small  compared  with  the  pressure  jump.  The  approximate 
theory  including  the  pressure  jump  predicts  the  general  behavior  of  the  data,  though 
of  course  the  data  do  not  correlate  in  terms  of  M  y  alone.  It  is  informative 

to  note  that  the  approximate  calculation  shows  that  there  should  be  significant  slip 


effects  (10%  or  greater)  in  the  weak  interaction  regime  M\/W;  **  10'“ and 
that  the  data  obtained  at  the  highest  densities  are  close  to  the  weak  interaction  limit. 

Discussion  of  the  Slip  Boundary  Conditions.  The  slip  boundary  conditions, 
given  by  Eqs.  14  and  17  are  usually  replaced  by  the  first-order  boundary  con¬ 
ditions;  that  is,  the  mean-free-path  is  evaluated  at  the  wall  conditions  instead  of 
conditions  at  the  base  of  the  boundary  layer.  The  validity  of  this  approximation 
can  be  tested  by  transforming  the  correct  boundary  conditions,  and  comparing  them 
with  the  first-order  boundary  condition.  With  the  following  transformation, 

ft.  ’ 0  ■  &&  ■  i  -  r 
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the  slip  boundary  conditions  are 


The  first-order  slip  boundary  conditions  obtain  if  one  replaces  — L_  with  unity. 

T* 

If,  however,  the  slip  boundary  conditions  are  again  substituted,  quadratic  expres¬ 


sions  in  and  ©Co)  are  obtained.  For  the  purpose  of  illustration, 

let  us  assume  local  similarity  so  that  (o)  =  O(o)  .  The  resulting  slip 


boundary  condition  is 


The  first-order  slip  boundary  condition  is  usually  obtained  by  noting  that  all  of 
the  terms  in  Eq.  (46)  are  of  higher  order  in  P*  than  the  last  term  under  the 


radical  so  that 


■  r>/¥ 


Examination  of  Eq.  (46)  in  light  of  the  cited  experimental  data  shows  that 

Tv 

this  boundary  condition  is  valid  only  for  the  case  of  an  adiabatic  wall,  ■■  — 1. 0. 

lo 

This  can  be  seen  by  invoking  the  strong  interaction  pressure  distribution  to  evaluate 
the  transformed  chordwise  variable. 
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For 


JL  =.  ‘K  Xe(n. +  i.'T37pi)'X  , 
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Now,  for  a  typical  cold-wall  case  St-sd 

*To 

The  experimental  data  show  that  the  low-density  departures  from  strong  inter- 
action  theory  begin  at  f”1  <a  0.  1,  and  all  of  the  terms  in  the  correct 

slip  boundary  condition,  Eq.  (46)..  are  of  the  same  order.  Consequently  it  is 
not  valid  to  consider  a  first-order  slip  boundary  condition  on  the  basis  of  ordering 
the  parameter  f”*  alone.  The  wall  temperature  effect  can  be  of  equal  impor¬ 
tance  and  the  entire  slip  boundary  condition  must  be  retained  for  the  case  of  a 
cold  wall. 

Comparison  with  Continuum  .Leading-Edge  Solutions.  In  the  discussion  on 
the  available  theories  it  was  noted  that  the  leading-edge  problem  had  been  treated 


by  Oguchi  with  continuum  concepts  assuming  slip  at  the  surface  could  be  neglected, 

and  taking  account  of  the  classical  (first-order)  velocity  and  temperature  jumps  at 
32 

the  surface.  The  basis  for  this  theory  was  originally  noted  in  Ref.  14  where  it 
was  shown  that  at  the  leading  edge,  the  flow  should  approach  a  wedge  flow;  that  is, 
a  straight  shock  wave.  This  model  is  treated  in  detail  in  Refs.  31  and  32  to  effect 
a  solution.  Note  that  in  those  developments  the  model  used  is  one  in  which  the 
shock  layer  is  fully  viscous.  However,  in  order  to  use  the  uniform-flow  Rankine- 
Hugoniot  relations,  it  is  required  that  the  shear  vanish  at  the  shock  wave.  This  is 
conceptually  difficult  to  accept  but  the  neglect  of  shock  wave  slip  does  not  appear 
critical  in  view  of  the  preceding  discussion.  In  addition,  the  pressure  jump  is 


31 
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neglected. 

The  slip  and  no-slip  solutions,  as  numerically  reported  in  Ref.  32,  are 
compared  with  the  present  experiments  in  Figs.  13  and  14  as  Ch  and 

j_  JL 

M2  Poo 

relation  is  not  obtained  within  these  parameters.  In  addition,  both  theories  given 
by  Oguchi  overestimate  the  heat  transfer.  This  could  stem  from  the  first-order 
boundary  conditions  which  underestimate  the  magnitude  of  the  slip  velocity  and 
enthalpy  jump  by  a  factor  of  2  to  3. 

An  interesting  item  to  be  noted  in  the  theory  of  Ref.  32  is  the  fact  that  a 

significant  slip  effect  on  heat  transfer  is  predicted.  This  is  in  contrast  to  other 

28 

results  for  zero-pressure  gradient  and  for  a  strong  interaction  pressure  grad- 
30 

ient  where  it  is  concluded  that  there  is  no  effect  of  slip  on  heat  transfer  because 
of  the  additional  sliding  friction  term.  The  fact  that  the  experiments  show  a  de¬ 
creased  heat  transfer  level  lends  partial  support  to  the  solutions  based  on  a  wedge 
flow. 

The  pressure  data,  Fig.  14,  show  the  same  general  features  as  the  heat 

transfer  data.  The  slip  and  no- slip  solutions  from  Ref.  32  appear  to  agree  with 

the  general  trend  of  the  data  obtained  at  the  higher  densities,  but  this  is  largely 

fortuitous  because  the  governing  parameters  do  not  produce  a  general  correlation. 

The  pressure  jump  given  by  Eq,  (24)  has  been  applied  to  the  slip  solution,  and  it 

can  be  seen  that  the  results  are  in  reasonable  agreement  with  the  data  obtained  at 

the  lowest  test  densities.  Again,  this  cannot  be  interpreted  as  a  confirmation  of 

the  theory  because  the  theory  does  not  produce  a  general  correlation. 

Comparison  with  Kinetic  Theory.  Solutions  have  been  obtained  by  Charwat, 

33  34,35 

using  approximate  kinetic  theory,  for  the  leading-edge  problem.  ’  ’  The 


Both  figures  show  again  that  a  universal  data  cor- 
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pressure  data  reported  here  are  not  rigorously  applicable  for  testing  this  theory 
since  the  measurements  were  made  5-100  ambient  mean-free-paths  from  the 
leading  edge.  They  do,  however,  provide  a  qualitative  check  on  the  theoretical 
trends.  The  heat  transfer  data  obtained  with  the  small-scale  model  provide 
detailed  information  throughout  the  transition  range  and  offer  a  quantitative  check 
on  the  theory. 

The  leading-edge  heat  transfer  data  are  presented  in  Fig.  15,  normalized  by 
the  theoretical  result  for  free -molecule  heat  transfer  (unit  accommodation  coef¬ 
ficient)  and  plotted  in  terms  of  distance  from  the  leading  edge  in  ambient  mean- 
free-path  length.  There  are  two  items  to  be  noted  in  these  data,  First,  both 
sets  of  data  indicate  that  in  the  limit  of  0  ,  the  thermal  accommodation 

coefficient  will  approach  unity.  Second,  the  heat  transfer  approaches  a  near¬ 
constant  value  several  mean-free-paths  from  the  leading  edge. 

One  point  of  potential  importance  in  interpreting  the  heat  transfer  data,  Fig. 

15,  is  the  fact  that  the  model  was  heated  by  the  current  in  the  thin-film  thermometers. 
Since  the  gages  were  not  uniformly  spaced,  the  model  had  a  nonuniform  axial 
temperature  distribution.  The  model  temperature  was  inferred  after  the  experi¬ 
ments  by  recording  the  gage  resistance  under  typical  pre-test  conditions.  These 
steady-state  surface  temperatures  are  presented  in  Fig.  16.  It  can  be  seen  that  the 
temperature  gradient  near  the  leading  edge  was  quite  severe  and  one  might  have  to 


In  this  connection,  it  should  be  noted  that  the  surface  was  titanium  dioxide  and 
no  unusual  attempts  were  made  to  clean  the  surface.  In  the  normal  course  of  the 
experiments,  the  model  was  maintained  at  a  temperature  of  about  150-250’F  and 
at  a  pressure  of  about  0.  1  -  1.0  pHg  for  one  to  two  hours  before  the  experiment. 
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take  account  of  thermal  creep  in  applying  these  data.  The  correction  due  to 
thermal  creep  has  not  been  included  in  any  correlations  presented  here,  It  has 
been  estimated  that  the  maximum  velocity  due  to  thermal  creep  is  about  5%  of 
the  free- stream  velocity. 

It  was  noted  in  the  discussion  of  the  theories  (Section  II)  that  the  approximate 

33 

kinetic  theory  for  the  leading  edge  problem  might  be  quantitatively  inaccurate 
owing  to  the  approximations  employed  in  accounting  for  the  incident  and  re-emitted 
molecules.  The  theory  should,  however,  identify  the  governing  parameters  and  pro¬ 
vide  a  basis  for  correlating  different  experiments.  The  heat  transfer  data  are  cor¬ 
related  in  terms  of  those  parameters  and  compared  witli  the  kinetic  theory  in  Fig. 

17.  It  can  be  seen  that  the  theoretical  parameters  correlate  the  data  from  the  two 
experiments  reasonably  well  with  a  scatter  of  about  8%.  In  the  free-molecule  limit, 
the  theory  overestimates  the  heat  transfer  by  a  factor  of  about  5.  If  the  theory  is 
forced  to  the  correct  free-molecule  limit,  as  in  Ref.  35,  it  is  found  that  the  theory 
predicts  the  heat  transfer  with  reasonable  accuracy  for  values  of  the  abscissa  of 
less  than  unity,  the  theoretical  range  of  validity. 

The  pressure  data  are  compared  with  the  kinetic  theory,  Fig.  18,  in  a  plot 


This  form  is  used  for  the  data  correlation  because 


of  _EL  vi 

Ro  •  i—  3  3 

is  the  basic  variable  in  the  kinetic  theory.  It  would  be  misleading  to  correlate 


the  pressure  data  in  a  manner  analogous  to  Fig.  17  because  the  pressure  data  are 
far  outside  the  range  of  validity  for  that  theory,  It  can  be  seen  in  Fig.  18  that  the 
present  data  are  reasonably  self-consistent  and  indicate  a  near-constant  pressure 


The  velocity  induced  by  the  higher  energy  molecules  emitted  from  the  hotter 
section  of  the  wall. 
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for 


»  10-100.  Included  in  Fig.  18  are  the  pressure  data  from  Ref.  36. 

It  can  be  seen  that  those  data  fall  some  50%-200%  above  the  present  data  and 
indicate  a  continuing  increase  in  pressure  in  the  range  ^  10-100.  The 
present  authors  have  been  unable  to  find  any  apparent  reason  for  this  discrepancy 
between  the  two  experiments  from  the  available  information. 


*- 
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CONCLUDING  REMARKS 

This  experimental  study  of  low  density  viscous  phenomena  at  near-free- 
molecule  conditions  has  shown  the  utility  and  versatility  of  both  the  hypersonic 
shock  tunnel  and  the  sharp  flat  plate  in  low  density  research.  Meaningful  measure¬ 
ments  in  millisecond  time  intervals  are  feasible  with  the  shock  tunnel  in  this  regime 
because  it  is  possible  to  use  high  reservoir  enthalpies,  and  hence  obtain  high  momen¬ 
tum  flux  and  energy  flux  in  the  test  section.  The  sharp  flat  plate  is  a  particularly 
useful  tool  because  many  of  the  competing  low  density  phenomena  are  suppressed 
with  this  elementary  configuration.  Hence,  the  governing  low  density  mechanisms 
can  be  isolated. 

The  experimental  data  presented  for  a  sharp  flat  plate  at  large  angles  of 
attack  are  in  good  agreement  with  viscous  shock  layer  theory  and  indicate  that  the 
small  discrepancies  between  theory  and  experiment  (10-15%)  can  be  attributed  to 
differing  values  of  the  Prandtl  number  in  the  theory  and  the  experiments.  Both 
theory  and  experiments  indicate  that  the  transport  effects  at  the  shock  wave  (shock 

I 

wave  slip)  are  negligibly  small  for  ~^gg  '  y  >°-  4- 

In  examining  the  applicable  theories  in  Section  II  and  V,  two  extensions  have 
been  made  here.  First,  there  is  a  pressure  jump  at  the  surface,  analogous  to  the 
classical  velocity  and  temperature  jump  at  the  surface.  This  effect,  which  is  a 
consequence  of  conserving  the  normal  component  of  momentum,  is  larger  than  the 
changes  in  induced  pressure  due  to  slip-induced  alterations  in  the  displacement 
thickness.  Evidently  this  effect  has  not  been  considered  in  boundary  layer  analyses 
that  include  the  slip  boundary  conditions.  The  second  extension  was  to  show  that 
the  first-order  slip  boundary  conditions  used  in  the  boundary  layer  analyses  are 
inapplicable  for  the  case  of  a  cold  wall.  In  writing  down  the  complete  boundary 
condition,  it  was  shown  that  since  the  wall  temperature  term  is  of  the  same  order 
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as  the  governing  parameter,  the  terms  usually  neglected  are  the  same  magnitude 


as  the  first-order  term.  Hence  for  ~  0.1,  the  first-order  slip 

To 

boundary  condition  is  a  factor  of  2  to  3  smaller  than  the  complete  boundary  con¬ 
dition. 


From  the  results  of  the  experiments  at  zero-angle  of  attack,  it  appears  that 
slip  at  the  surface  is  the  dominant  .  luse  of  the  observed  low  density  departures 
from  the  thin  boundary  layer  regime.  This  conclusion  is  consistent  with  the 
observations  on  the  pressure  jump  and  the  inapplicability  of  the  first-order  boun¬ 
dary  conditions.  The  pressure  jump,  when  applied  to  existing  theory,  brings  theory 
and  experiment  into  qualitative  agreement.  The  discrepancies  between  experiment 
and  theory  for  heat  transfer  might  well  be  due  to  the  fact  that  the  first-order  boun¬ 
dary  conditions  underestimate  the  velocity  jump  and  enthalpy  jump.  The  other 
low  density  effects  are  undoubtedly  present;  however,  their  combined  effects  are 
evidently  small  in  comparison  with  that  due  to  surface  slip.  Further  studies  are 
required  to  quantitatively  assess  the  other  low  density  effects. 

The  present  heat  transfer  results  define  the  transition  to  the  free-molecule 
limit,  but  further  experiments  over  a  wider  range  of  conditions  and  closer  to  the 
leading  edge  would  be  of  value.  The  most  pressing  experimental  need,  however, 
is  for  high  Mach  number  pressure  data  close  to  the  leading  edge.  The  available 
pressure  data  are  not  adequate  to  test  the  kinetic  treatments  of  the  leading  edge 
problem. 
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Figure  2  SMALL-SCALE  FLAT  PLATE  MODEL 


Figure  4  LARGE-SCALE  FLAT  PLATE  MODEL 


Figure  5  HEAT  TRANSFER  OH  A  SHARP  WEDGE  TQ  »  3600 “K 
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Figure  6  TRAIUHG  EDGE  EFFECTS  III  WEDGE  FLOW 
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Figure  8  PRESSURE  D I STR I  BUT I  ON  ON  A  SHARP  FLAT  PLATE 


Figure  II  PRESSURE  DISTRIBUTION  ON  A  SHARP  FLAT  PLATE 


Figure  15  HEAT  TRANSFER  HEAR  THE  LEADING  EDGE  OF  A 
SHARP  FLAT  PLATE 


